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We show that the triangular-lattice antiferromagnet Q-NaMn02 is endowed with a complex mag- 
netostructurally inhomogeneous ground state, a result of intertwined spin and lattice degrees of 
freedom. Synchrotron X-ray diffraction, nuclear magnetic resonance and muon spin relaxation un- 
veil that the magnetically frustrated monoclinic phase is disrupted by nanoscale regions of a triclinic 
structure. The magnetic order in the latter is uncorrelated with the monoclinic matrix, where bulk 
Neel order is initiated by the local-scale inhomogeneity. The phase separation is endorsed by frus- 
tration leading to quasi-elastic magnetic excitations that couple to competing structural phases. 

PACS numbers: 75.80.+q, 76.60.-k, 76.75.+i, 64.75.-g 
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Complexity of phases is a common theme in systems 
with simultaneously active degrees of freedom - spin, 
charge, orbital, and/or phonon [7|. Owing to strong elec- 
tronic correlations, nanometer-scale inhomogeneities, like 
stripes and checkerboard charge ordering in cuprates 
01 , as well as maiiy other striking inhomogeneous states 
of manganites [a , cobaltates nickelates d, 0|, and 
iron pnictides |10l lll| , can emerge even when Hamiltoni- 
ans do not brake translational symmetry 0, [l2|- In these 
systems inhomogeneity stems from the charge degrees of 
freedom and usually appears due to quenched disorder 
- random deviations from perfect system uniformity - 
that can dramatically alter the properties of materials en- 
tailing competing interactions [7|. In localized-spin sys- 
tems, on the other hand, inhomogeneous phases are rare 
0, [3| and still poorly understood, although they have 
been predicted to appear already in the clean limit of 
uniformly frustrated ma gnet s [T El [l6| . Geometrical frus- 
tration of a spin lattice [17|, ll8| is known for provoking 
novel fundamental phenomena in magnetism, like exotic 
disordered ground states (isl . [l9| and fractional spin exci- 



tations 



Moreover, its involvement in structural 
instabilities has also been well documented, e.g., in re- 
laxor ferroelectrics [2l| and multiferroics 22 1. Therefore, 
ground-state degeneracy of geometrically frustrated mag- 
nets, in alliance with the magnetoelastic coupling, could 
make the energetic cost of spatial fluctuations low, ren- 
dering these systems "electronically soft" [S^ and thus 
prone to forming inhomogeneous states. 

The layered a-NaMn02 compound is intriguing in this 
respect, because of its atypical low-temperature magnetic 
state. It represents a spatially anisotropic triangular lat- 
tice of S' = 2 (Mn'^+) spins with dominant intrachain 
(Ji = 65 K) and frustrated interchain (J2 = 29 K) an- 
tiferromagnetic exchange interactions [2| in the ah crys- 



tallographic plane (Fig. [T]). Both, sharp magnetic re- 
flections and diffuse scattering were simultaneously ob- 
served in neutron powder diffraction (NPD), even far 
below the Neel temperature Tjv = 45 K [J. The co- 
existence of long-range magnetic order and short-range 
spin correlations indicates a magnetically inhomogeneous 
ground state. In addition, large and highly anisotropic 
microstrains in the crystal structure above Tjv suggest 
structural inhomogeneities. This is a likely precursor to 
the structural phase transition from a monoclinic to a 
triclinic symmetry, which was reported to occur concur- 
rently with the magnetic transition due to magnetoelastic 
coupling [H . Whether structural inhomogeneities remain 
present below T^v and if they pertain to a possible mag- 
netic phase separation are major unresolved issues. 




FIG. 1: Anisotropic triangular lattice of S = 2 (Mn^^) spins 
in a-NaMn02 with exchange interactions Ji, J2 in the ab 
crystallographic plane of the monoclinic (m) and the triclinic 
it) cell settings. Long-range antiferromagnetic order [km = 
(^,^,0)] is interrupted by the shaded nano-sized magnetic- 
defect region, characterized by a correlation length ^. 
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In this Letter, we combine synchrotron X-ray powder 
diffraction (XRD) with local-probe nuclear magnetic res- 
onance (NMR) and muon spin relaxation (/iSR) tech- 
niques, which provide insight to a unique magnetostruc- 
turally inhomogeneous ground state of a-NaMn02. Its 
spatial complexity stems from simultaneously active spin 
and lattice degrees of freedom. We show that the re- 
lief of frustration by a lattice-distorting procedure, likely 
triggered by infinitesimal quenched disorder below Tjv , is 
endorsed by quasi-elastic magnetic excitations and occurs 
with the spontaneous emergence of competing magnetic 
states on nanoscopic level. 

We have performed a thorough structural investi- 
gation of a-NaMn02, employing high- resolution syn- 
chrotron powder XRD (ID31 beam-line at the ESRF; 
A = 0.39986 A, T = 5-300 K) to gain insight on the large 
strain fields detected in the early NPD study [ij. The 
XRD confirms the single-phase nature of a-NaMn02. 
At room temperature it crystallizes in the monoclinic 
structure [C2/m; a = 5.67075(5) A, 6 = 2.85896(2) A, 
c = 5.80109(6) A, /3 = 113.1537(4)°; Fig.pta)]. Rietveld 
refinements with the FULLPROF suite 26| show that, 
upon cooling, /ifcZ-dependent anisotropic broadening of 
the Bragg refiections needs to be taken into account. We 
have utilized the Stephens's formalism Q of anisotropic 
strain broadening (Shki) and found that the iS'o4o param- 
eter is the largest and increases the most down to Tpf 
[Fig. [5I^b)]. This suggests that the distribution of the ac- 
interplane distances is the broadest and highly affected 
when approaching Tjv from above, while the ab (6*004) 
and be (6*400) interplane-distance distributions are much 
smaller and less temperature dependent. Qualitatively, 
the iS'o4o strains in the monoclinic structure may origi- 
nate from minute shifts of the Mn-Mn spin-chains along 
the bm axis (see Fig. [T] and Ref. 1281) . a lattice distortion 
which can effectively lead to the triclinic structure (PI) 
at T < T/v Q . This behavior implies that the system de- 
velops a precursor state, before the magnetoelastic effects 
lift the degeneracy of its ground state. 

In contrast to the isostructural CuMn02 @|i the 
high-resolution synchrotron XRD patterns of Q;-NaMn02 
below Tpf do not show any resolvable peak-splitting that 
would unambiguously confirm a bulk phase transforma- 
tion to the lower-symmetry triclinic phase. Surprisingly, 
Rietveld refinements at T < T/v on the highly parame- 
terized, strained triclinic (t) phase produce a quality of 
fit which is statistically poorer compared to the model 
based on the strained monoclinic (to) phase, Xt ~ l-^Xm 
[ist . Furthermore, a close inspection of the Rietveld 
fitted patterns for the latter recognizes anisotropically 
broadened profiles, where the calculated pattern system- 
atically misses intensity from the low- and high-angle 
tails of certain Bragg reflections [Fig. [2fc)]. It is worth 
noting that presence of defects [3l|j chemical or struc- 
tural inhomogeneities (3^ may be the reason for the de- 
velopment of short-scale regions, which can give-rise to 
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FIG. 2: (a) Observed (circles), calculated (full curve) and dif- 
ference synchrotron XRD profiles of a-NaMn02 at 295 K, af- 
ter Rietveld refinement in the monoclinic symmetry (C2/m); 
predicted reflections are marked by vertical ticks. At 26 > 21° 
all pattern intensities were multiplied by 3 to increase clar- 
ity, (b) The temperature evolution of the anisotropic strain 
broadening parameters Shki (lines are a guide to the eye). 
Detail of the 30 K proflle after Rietveld flts with (c) a single 
strained monoclinic phase (important deviations in the differ- 
ence plot are marked by arrows) and (d) coexisting strained 
monoclinic and unstrained triclinic phases. 



unusually broadened profiles, like those of a-NaMn02. 
A way of simulating such a behavior is by a two-phase 
(2p) Rietveld model accounting for a nanoscale mixture 
of phases. Significantly improved quality of fit [Fig.[5Jd)[, 
xip ~ 0-35Xmi is obtained with a spatially inhomoge- 
neous model assuming coexistence of the monoclinic and 
the triclinic phases [281 . The latter has been found to 
attain a volume fraction of ~35(1)% at 30 K, which 
marginally evolves down to 5 K. The majority mono- 
clinic phase retains some strain, but this is severely re- 
lieved below T/v [Fig. [l{b)[ where the triclinic phase, 
with inequivalent interchain Mn-Mn distances |28| . co- 
exists with the former. Randomly distributed nanoscale 
inhomogeneities, pertaining to the minority phase (PI) 
that is grown at the expense of the majority one (C2/to), 
may severely modify the macroscopic electronic response 
(vide infra) of the bulk compound. This renders it an 
intriguing subject for local-probe investigations. 

In order to validate the apparent phase separation and 
associated local inhomogeneities in a-NaMn02, we em- 
ployed 23Na (/ = 3/2) NMR [l^, locally probing elec- 
tronic (charge and magnetic) environment of the ^'^Na 
nuclei. A transferred hyperfine coupling A,/, of the nu- 



clear spin Ij with electronic spins Sk on nearby Mn'^+ 
sites, ~ '^j'J2k—jk ' ^A; = /i7lj-5-Bo (fi is the reduced 
Planck constant and 7 = 70.8 MHz/T the ^^Na gyromag- 
netic ratio) . is responsible [2^ for a sizable isotropic shift 
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FIG. 3: (a) The central (-i o i) line of the ^^Na NMR 
spectrum in a-NaMn02 at 300 K and 8.9 T. The solid line 
corresponds to the fitting, explained in the text, (b) Scaling 
of the isotropic Siao (line shift) and anisotropic 5aniso (line 
width) hyperfine coupling constants with bulk susceptibility 
X at 5 T. Solid lines are linear fits of the high-temperature 
data, (c) Partial wipe-out of the NMR signal (corrected for 
spin-spin relaxation) below Tn- (d) The ^'^Na spin-lattice 
relaxation rate exhibiting power-law dependence below Tn- 
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FIG. 4: (a) The decay of /i^ polarization above and below 
Tat = 45 K. The solid line corresponds to exponential decay, 
(b) The temperature evolution of the local magnetic field at 
muon sites (lines are a guide to the eye), (c) Experimental os- 
cillations of /x^ polarization in the ordered state (circles). The 
lines show the fitting with two damped cosine contributions 
and modeling with nanoscopic inhomogeneities, as explained 
in the text. Inset: contour of the electrostatic potential in the 
ab plane with highlighted global minima for the Pi site. 



[6iso = 707(10) ppm at 300 K] and anisotropic broadening 
[^aniso = 950(10) ppm] of the central (— ^ O i) nuclear 
transition [Fig. EI a)]. Our calculations [2^ show that 
<5aniso corrcsponds to local dipolar fields while the finite 
i5iso proves that the inter-layer exchange effectively runs 
through the Na+ sites. 

The (5iso and (5aniso are determined by the average of 
a static local-field distribution at the ^'^Na site and its 
width, respectively. Therefore, they are both expected to 
scale with magnetic susceptibility x- This should remain 
the case even below T/v, since on the highly symmetric 
^■^Na position (structural center of inversion) the local 
fields from antiferromagnetically ordered Mn"^"*" moments 
at sites ±r [J effectively cancel out. Experimentally, 
we observe significant deviations form the linear scal- 
ing, particularly pronounced for i^aniso [Fig. E^b)], which 
set in already at ~ 60 K, i.e., far above Tn — 45 K. 
Since the ^'^Na quadrupolar-broadening effects are small 
[2^ the deviations above Tn signal an onset of short- 
range spin correlations. On the other hand, the apparent 
non-linearity (Janiso 9^ X below T/v discloses a more com- 
plicated magnetic order than initially proposed [H. We 
stress that (5aniso does not exhibit an order-parameter de- 
pendence below T/v, hence, it cannot be explained by a 
small additional spin component that would break the 
symmetry at the ^■^Na sites. Quite importantly, we also 
find that close to Tn the NMR intensity of the cen- 
tral line that would normally scale with 1/T due to the 
Boltzmann population factor, suddenly drops by 30% 



[Fig-ISIc)]- NMR thus clearly discloses a magnetic phase 
separation below T/v , where the remaining NMR compo- 
nent characterized by a zero static local magnetic field is 
accompanied by an unobservable component. 

In order to expose the microscopic nature of the mag- 
netic phase separation below T/v, we finally resort to 
/iSR. This technique is extremely sensitive to local mag- 
netic fields B^, and directly measures their distribu- 
tion in an ordered state through the relaxation func- 
tion P{t) of the muon polarization 0. Zero-field (ZF) 
/iSR experiments were performed on the General Pur- 
pose Surface-Muon Instrument (GPS) at the Paul Scher- 
rer Institute, Switzerland. The phase transition at Tn 
is clearly visible in reduced muon polarization below T/v 
[Fig. |ll^a)[. In a static local field such decrease re- 
sults from averaging out fast oscillations in the relax- 
ation function J§[ T'pwd(i) = 5 + |cos(7^i3pt) of a pow- 
der sample. We stress that below T/v the polarization 
drops to 1/3, which is a clear indication that the bulk 
of the sample has undergone a transition into the or- 
dered state. Oscillations in P{t) below T/v are evident 
at early times [Fig. |lllc)[. They reveal that two in- 
equivalent muon stopping sites exist, as the polariza- 
tion can be nicely fitted to the two-component model 
Pit) = EU [i^"^'* + |cos(7Mi?Mj)e-^*-*]. We find 
volume fractions /i:/2 = 60(5):40(5) that are temper- 
ature independent and internal fields j that show an 
order-parameter-like temperature dependence [Fig.|31[b)]. 

We emphasize the large difference between the Ion- 
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gitudinal A; = 0.03(1) /is~^ and the transverse Ati = 
82(5) ns^^; At2 ~ 12(1) fis^^ relaxation rates at 5 K. 
The former is of dynamical origin, while the latter 
additionally encompasses a distribution of static local 
fields Q. This distribution is obviously responsible for 
the fast damping of oscillations. In order to model 
it, the knowledge of muon stopping sites is required. 
Here, we rely on ab-initio calculation of the electro- 
static potential of a-NaMn02 [l^, which predicts that 
the muons most likely stop at the 4z interstitial sites 
Pi = (0.18,0,0.47) and P2 = (0.45,0,0.75), the former 
corresponding to the global electrostatic-potential mini- 
mum [inset in Fig. [41(c)]. The calculated dipolar fields 
from antiferromagnetically ordered Mn'^"'' moments on 
these sites match the experimentally determined fields 

[Hi. 

The broad local-magnetic-field distributions at Pi and 
P2 provide a microscopic insight to the inhomogeneous 
magnetic state of Q;-NaMn02. Most importantly, they 
can be related to the "parasitic" triclinic structural 
phase. We model the field distributions by introducing 
triclinic magnetic defects in the form of nano-sized is- 
lands with the magnetic order reversed at the interfaces 
with the monoclinic matrix (Fig. [1]). We have consid- 
ered these defect regions to be either one- (ID), two- 
(2D), or three-dimensional objects with a characteristic 
correlation length ^, randomly positioned in the lattice 
so that a given fraction p of Mn'^+ ions resides in them. 
Independent on their dimensionality, when ^ = 00 (corre- 
sponding to stripe- and slab-shifts for ID and 2D defects, 
respectively) we find narrow local-field distributions with 
only a few characteristic frequencies [2^. These are man- 
ifested in beating [Fig. Hl^c)] rather than in damping of 
P{t). On the other hand, for finite ^ the field distri- 
butions become much broader [2^ and, consequently, 
result in heavily damped P{t). For p ~ 0.35, as de- 
rived from the above-presented structural analysis, we 
can reproduce the experimental damping of P(i) at 5 K if 
^ = 30(5) A [Fig. m^c)]. Such dispersed finite-size defects 
are in conformity with the short-range ordered regions in- 
ferred by neutron diffraction P). Moreover, introduction 
of such defects also explains the sudden ~30% decrease 
of the NMR intensity below Tjv. Namely, the defects 
locally break the inversion symmetry and lead to a two- 
peak local-field distribution at the ^"^Na sites [2^. The 
dominant peak remains positioned at BNai = and is 
modestly broadened, which may explain the increase of 
(^aniso bclow ~ 60 K [Fig. [3]^b)] . The second peak, on the 
other hand, is found at i3Na2 = 0.26 T, yielding the huge 
width of the corresponding NMR signal i?Na2/Bo = 3% 
that renders this component unobservable. 

In order to shed light on the mechanism responsible for 
stabilizing the remarkable inhomogeneous ground state 
of Q!-NaMn02, we last inspect its low-energy magnetic 
excitations. According to inelastic neutron scattering, a 
large gap (Aq — 87 K) due to single-ion anisotropy and 



ID excitations dictate the ground state at zero field 
Since these measurements can easily miss a quasi-elastic 
part, we resort to the complementary ^^Na spin-lattice 
relaxation study. Below T/v, the relaxation is governed 
by inelastic scattering of spin waves [3^. For gapped ex- 
citations the thermally-activated relaxation rate 1/Ti oc 
g-A/T jg expected for T <C A, and the power-law depen- 
dence l/Ti oc T" for T > A [Hi. Here n > 2 depends 
on the dimensionality of the spin-wave dispersion and 
the number of magnons assisting the relaxation process 
[ii,!!!]. The sharp suppression of l/Ti at T/v [Fig.[3)^d)] 
proves that NMR indeed probes low-energy magnetic ex- 
citations. However, the expected exponential law with 
A = Aq — /iBq = 53 K (/X = 2.9^B \M) is not in consent 
with the data [Fig. [3jd)]. In clear contrast, the power 
law with n ~ 4.3 (and the additional low-temperature 
constant term 1/rf = 0.01 s~'^ corresponding to resid- 
ual relaxation) matches the experimental data very well. 
Deviations from this law are observed only in a critical re- 
gion close to Tat. Thus, our NMR spin-lattice relaxation 
measurements clearly witness quasi-elastic magnetic ex- 
citations. These can be explained by sliding magnetic 
stripes along the bm axis, an effect that costs zero energy 
in the monoclinic setting due to the frustrated nature of 
the interchain exchange J2 . As the individual stripes can 
shift independently, the ensemble of such excitations be- 
comes extended. In alliance with the magnetoelastic cou- 
pling, they provide a likely rationalization of the observed 
unique magnetostructural inhomogeneity in a-NaMn02. 

In summary, we have discovered a magnetostructurally 
inhomogeneous ground state in a-NaMn02 that entails 
defect nanoscale regions of the triclinic phase, where mag- 
netic order is altered with respect to the long-range or- 
dered monoclinic phase. The Neel order in the latter is 
likely governed by interphase boundaries or structural 
defects which pin the spin arrangement in the other- 
wise frustrated triangular layer topology. The observed 
inhomogeneity is in sharp contrast to superconducting 
cuprates and iron pnictides, colossal-magnetoresistant 
manganites and the majority of other examples, where 
microscopic phase separation originates from charge de- 
grees of freedom [2t3, [2"12|. Moreover, the existence 
of nanodomains in these systems depends on doping, 
while a-NaMn02 resides in the clean limit. We asso- 
ciate its remarkable ground state to a fierce competition 
between the monoclinic and the triclinic structure, whose 
energies are indistinguishable according to ab-initio cal- 
culations [33. The energy difference remains minimal 
when exchange interactions are considered 37, 3^. This 
near-degeneracy enables even infinitesimal quenched dis- 
order to locally lift the inherent geometrical frustration of 
the triangular lattice in the monoclinic phase. Thus, we 
propose that the geometrical frustration can lead to inho- 
mogeneities in the magnetoelastic channel, which repre- 
sents an intriguing novelty among phase-separated states 
of highly correlated electron systems. 
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Structural Considerations 

High-resolution synchrotron XRD patterns were col- 
lected at the ID31 beam-line in the ESRF (A = 
0.39986 A, T = 5 - 300 K), on the same high-quality 
polycrystalline samples that were used before [l|-|3|. At 
300 K a-NaMn02 crystallizes in the monoclinic structure 
[C2/m; Rietveld refinement R- factors not corrected for 
background, Rp = 11.6%, i?wp = 15.8%, i?cxp = 3.80%, 
= 17.3; see Fig. 2(a) of the main part], a variant of 
the prototype rhombohedral delafossite, due to the pres- 
ence of Jalm- Teller active Mn^+ cations. Upon cooling, 
this phase evolves by showing /ifcZ-dependent anisotropic 
broadening, which in the Rietveld refinements was mod- 
eled by the Stephens's formalism {Shki) including 
a pseudo-Voigt peak-shape function. We find that the 
S'o4o parameter, measuring the covariance {{a2 — (0^2))^) 
of the distribution of the metric parameter a2 = 1/6*^ 
is the largest and increases the most between 300 K 
and Tjv [Fig. 2(b) of the main part]. The strain fields in 
the monoclinic structure can originate from minute shifts 
of the Mn spin-chains along the bm axis (Fig. [5]) . 

Although the ID31 synchrotron experiment provides 
an increased angular resolution (nominal instrumental 
contribution to the FWHM is 0.003° in 29) as com- 
pared to the previous NPD structural study (FWHM 
was ^ 0.3° in 26), the synchrotron XRD patterns 




FIG. 5: Schematic of the lattice-distorting in a-NaMn02 at 
the transition from the monoclinic (m) to the triclinic (t) 
phase. The Mn-Mn antiferromagnetic 5 = 2 chains shift along 
the bm axis. Arrows indicate the ordered magnetic moments. 



of a-NaMn02 below Tpf do not show any resolvable 
peak-splitting (Fig. [6]), in contrast to the isostructural 
CuMn02 phase [5|, |6|. In the latter compound such a 
feature signals a bulk phase transformation to the tri- 
clinic cell. Rietveld refinements at T < Tjv on the 
highly parameterized (as defined by the much larger num- 
ber of generated reflections due to -1 Lane symmetry) 
strained PI phase produces a quality of fit [Rp — 17.2%, 
i?wp = 22.4%, i?exp = 3.9%, = 32.3; Fig. [TJa)] 
which is poorer than that of the strained C2/m phase 
[Rp = 14.9%, i?wp = 18.9%, i?cxp = 3.9%, ^ 23.1; 
Fig. Elb)]. The a-NaMn02 data ai T < Tn lack peak- 
splitting, however, a close inspection of the Rietveld fit- 
ted patterns in the monoclinic phase recognizes that the 
calculated profile systematically misses some intensity 
from different regions of the observed diffraction profile 
[Fig. Ha)]. 

Modeling such anisotropically broadened profiles may 
assume a nanoscale mixture of phases, where, in a sim- 




8.0 8.5 9.0 9.5 8.0 8.5 9.0 9.5 



2e{°) 2e(°) 

FIG. 6: Part of the synchrotron powder X-ray diffraction pat- 
terns from the ID31 beam-line at the ESRF (A = 0.39986 A) 
for (a, b) a-NaMn02 and (c, d) the isomorphous CuMn02, 
above and below the transition to the Neel ordered state (Tjv). 
The /ifeZ-dependent anisotropic broadening of the Bragg peaks 
is evident in the former, whereas characteristic peak-splitting 
of certain families of Bragg reflections in CuMn02 signifies a 
bulk phase transformation from the monoclinic to the triclinic 
structure. 
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FIG. 8: A detail of the 30 K profile of Q-NaMn02 fitted by 
the Rietveld method on the basis of (a) a single-phase model 
entailing a strained monoclinic phase and (b) a two-phase 
model involving the coexistence of strained monoclinic and 
unstrained triclinic phases. Arrows mark important devia- 
tions in the difference plot due to inadequacies of the calcu- 
lated profile of the single-phase model. 



FIG. 7: Single-phase Rietveld refinement of a-NaMn02 at 
30 K, (a) in the triclinic cell setting (a = 3.f637f(6) A, 
b = 3.f6f27(8) A, c = 5.775f9(6) A, a = 110.340(2)°, 
P = 110.316(1)°, 7 = 53.7094(8)° and (b) in the mono- 
clinic ceU setting, a = 5.64258(7) A, b = 2.85725(3) A, 
c = 5.77284(9) A, /3 = 112.9123(7)°. Large deviations be- 
tween observed and calculated profiles are depicted in the 
difference plot. 



plified approach, a multi-phase Rietveld fit describes the 
data. Once the strained monoclinic and unstrained tri- 
clinic phases are allowed to coexist, a significantly im- 
proved quality of fit is attained {Rp = 8.7%, i?wp = 
11.3%, i?exp = 3.9%, = 8.2; Fig.Hb), Fig.©, yield- 
ing a "parasitic" triclinic phase with a volume fraction 
of ^ 35(1)%, which is retained down to 5 K. Nanoscale 
structural inhomogeneities may unveil the presence of 
quenched disorder However, this cannot be well- 

resolved from the present powder diffraction data and re- 
quests high-resolution transmission electron microscopy, 
as exploited before for the non-stoichiometric CuMn02 
phases j^. The temperature evolution of the in-plane 
Mn-Mn separation indicates a lattice distortion proce- 
dure, which gives rise to the relief of strain and the in- 
equality of the two interchain Mn-Mn distances below 
about Tjv [Fig. [TOT a)] . The temperature evolution of the 
a-NaMn02 unit-cell volume shows normal thermal con- 
traction, whereas the coexisting triclinic phase holds half 
of the monoclinic cell volume [Fig. [TUTb)]. 



^^Na NMR experiments 

Frequency-swept ^^Na (J = 3/2) NMR spectra were 
measured on a home-built spectrometer in the fixed 
magnetic field Bq ~ 8.9 T and obtained with a solid- 
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FIG. 9: Observed (circles), calculated (full curve) and dif- 
ference synchrotron XRD profiles of Q-NaMn02 at 30 K, 
after Rietveld refinement with a two-phase model assuming 
coexistence of monoclinic (C2/m; a = 5.63750(3) A, b = 
2.85763(1) A, c = 5.77295(3) A, /3 = 112.8926(4)°) and tri- 
clinic (PI; a = 3.1814(2) A, 6 = 3.1751(2) A, c = 5.7589(2) A, 
a = 110.429(3)°, /3 = 110.589(4)°, 7 = 53.308(3)°) phases; 
predicted refiection positions are marked by the vertical ticks 
(top set : monoclinic and bottom set: triclinic space groups). 
At 28 > 21° all pattern intensities were multiplied by 3 to 
increase clarity. 
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FIG. 10: (a) Temperature evolution of the inplane Mn-Mn 
distances in the two unit- cell settings, monoclinic ((72/m) and 
trichnic (PI), as derived by the Rietveld refinements of the 
synchrotron XRD data. Below Tjv Rietveld fits were under- 
taken with the model assuming coexistence of strained mon- 
oclinic and unstrained triclinic phases, (b) The T-evolution 
of the monoclinic and triclinic unit cell volumes. 



echo pulse sequence. The 7r/2 pulse length was 3.4 /is. 
Spin- lattice relaxation (1/Ti) measurements utilized an 
inversion-recovery sequence at the peak of the central 
(— i o i) nuclear transition. The NMR line shift 
was measured relative to the reference Larmor frequency 
i/Q = 100.5255 MHz, corresponding to the NMR signal of 
a 0.1 M NaCl solution. 

The ^'^Na spectrum recorded at 300 K is shown in 
Fig. [TT] It consists of a sharp central line (—1/2 <-!> 1/2) 
and satellite transitions (±1/2 -o- ±3/2), broadened by a 
quadrupolar coupling in first order. The satellites and the 
central transition are strongly non-symmetric, implying 
that anisotropic magnetic coupling with Mn'^^ magnetic 
moments is large. The model Hamiltonian used for fitting 
the spectrum included the Zeeman part, the quadrupo- 
lar coupling of the nuclear quadrupolar moment with lo- 
cal electric-field gradients (EFG) and uniaxial anisotropic 
hyperfine coupling of the nuclear magnetic moments Ij 
and Mn'^"'' electronic moments, = hj J^j Ij '^'Bo {h is 
the reduced Planck constant, 7 = 70.8 MHz/T). Isotropic 
powder averaging was performed. The obtained EFG 
asymmetry parameter is 77 = 0.54(3) and the quadrupo- 
lar frequency is vq ~ 1.51(5) MHz, which is a rather 
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FIG. 11: (a) The ^^Na NMR spectrum of Q:-NaMn02 mea- 
sured at 300 K in the applied field of 8.9 T {vq = 
100.5255 MHz), (b) The central transition (-1/2 o 1/2) of 
the same spectrum. Solid lines correspond to the fit explained 
in the text. 



typical value for the ^^Na nuclei. The isotropic part of 
the hyperfine magnetic coupling tensor 5;^^, which is re- 
sponsible for the shift of the central line from the ref- 
erence frequency 1^0, is Siso = 707(10) ppm at 300 K. 
The anisotropic part responsible for the line shape of the 
central transition, has the principal eigenvalue (5aniso ~ 
950(10) ppm at 300 K. The quadrupolar interaction leads 
to the second-order shift -20 ppm and the width 100 ppm 
of the central line, which is small compared to the effect 
of the hyperfine interaction. The typical chemical-shift 
range -60 - 10 ppm for ^^Na is also negligible. 

Temperature dependence of i5iso and (5aniso [Fig. I12| 
provides a local insight to magnetism in a-NaMn02. 
Their expected scaling with magnetic susceptibility, 
<5iso, ^aniso oc X, is obscrvcd at high temperatures [cf. 
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FIG. 12: The temperature dependence of the ^^Na NMR line 
shift (5iso and the line width 5aniao, compared to bulk suscep- 
tibility X (solid line) measured at 5 T. 
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Fig. 3(b) of the main part]. The scahng of 6iso is roughly 
retained down to 15 K where the low-temperature up- 
turn occurs in X [H • The intrinsic susceptibility Xi cx ^iso 
continues decreasing down to the lowest temperatures. 
Much larger deviations of Janiso from linearity are ob- 
served already at ~ 60 K and suggest a more complicated 
magnetic state than initially suggested 

Magnetization-recovery curves in spin-lattice re- 
laxation measurements were fitted to Mz{t) = 
A{1- iexp [-(t/T^i)"] - fexp [-(6t/Ti)"]} + B, suit- 
able for magnetic relaxation of the central line of a 
/ = 3/2 nuclei Q. We found constant stretch expo- 
nent a = 0.86(2) at T > Tm and a significantly de- 
creased exponent in the magnetically ordered state reach- 
ing a = 0.60(2) at low temperatures. The decrease of 
a below Tn evidences an increased distribution of spin- 
lattice relaxation rates, hence a broad distribution of lo- 
cal environments, in line with the inhomogeneous ground 
state of a-NaMn02. 



Muon Stopping Sites 

It is supposed that would most likely stop at 
the global potential minimum of the crystal structure 
911, which has been verified before on several instances 
id . The calculations were carried out by apply- 

ing the Quantum Espresso code [l2|, and b y u sing the 



generalized-gradient approximation (GGA) [13| for the 
exchange-correlation potential. The electron-ion interac- 
tions were described with the Vanderbilt ultrasoft poten- 
tials The plane- wave cut-off parameters were set to 



408 eV and 3264 eV, respectively, whereas a 4 x 4 x 4 
(4 X 8 X 4 for the non-spin-polarized case) mesh in the 
full Brillouin zone was used for the Methfessel-Paxton 
sampling [lH integration. Both the spin-polarized and 
non-polarized calculations yield the same most-probable 
muon site Pi = (0.18,0,0.47) [inset in Fig.4(c) of the 
main text], which is an interstitial site with site symme- 
try m at Wyckoff position Ai in the monoclinic setting. 
Although this site is close to the Mn-layer, its nearest- 
neighbor site is actually at the distance of 1.74 A. 

Calculations of the dipolar field at this site {vide in- 
fra) yield the magnetic field of 0.24 T that nicely corrob- 
orates the experimental B^i = 0.22(2) T. However, our 
ab-initio calculation could not find any local minimum 
that would correspond to = 0.69(2) T. Therefore, 
our assignment of the second site relies on a common 
tendency of the muon to "bind" to at the distance 
of 1.0 A 0. We then set the second muon stopping site to 
P2 = (0.45,0,0.75) (4i) on the 1.0 A-sphere around the 
oxygen site, where the dipolar field matches i?p2 and the 
electrostatic potential of the unperturbed crystal struc- 
ture is minimal. 

The influence of a moun stopping at either Pi or P2 
on all of the exchange constants and thus locally on the 



magnetic ground state is expected to be negligible. The 
site Pi is far from any exchange path, while the site P2 
could possibly interfere only with the weak inter-layer 
path. Therefore, it is highly unlikely for the muon to 
cause any local inhomogeneities of the magnetic order. 



Dipolar-field calculations 

Calculations of the dipolar magnetic field Bdd(r) = 

BL V . - 3(r-r )[p,.(r-r,)] | ■ ■ 

Ait ^—ij I |r— rj|-* |r— rj|-> J " 

performed by taking into account all Mn'^''' moments ^.j 
at sites r^, within a sphere centered at r with a ra- 
dius large enough to ensure convergence. The moments 
were either taken polarized along the applied field in 
the paramagnetic state or ordered according to Ref. [l| 
(/i = 2.9 /is at 5 K). We find that the dipolar magnetic 
field at the ^^Na site in the paramagnetic phase is al- 
most axially symmetric. Its principal value amounts to 
A = 0.16 T/hb- Taking into account the 300 K suscep- 
tibility value X = 3.5 • 10~^ emu/mol, this leads to the 
dipolar contribution to the anisotropic hyperfine coupling 



Ax 
Na 



1000 ppm {Na is the Avogadro number). 
This matches astonishingly well with the experimental 
(^aniso — 950 ppm, which disclosed a negligible contribu- 
tion of the transferred hyperfine interaction to (5aniso- 

When modeling the normalized local magnetic-field 
distributions D{B) = at both muon stopping sites 

and at the ^^^Na site in the structurally inhomogeneous 
phase below Tat, we assumed that the magnetic order in 
the triclinic "parasitic" regions was inverted with respect 
to the monoclinic matrix. We calculated the magnetic 
field in No points, each with a different random distribu- 
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FIG. 13: Field distributions in the inhomogeneous magnetic 
phase of a-NaMn02, calculated for the defect concentration 
p = 0.35 on both muon stopping sites. Pi and P2, and on the 
^■^Na nuclear site. Results are shown for the one-dimensional- 
defect model with ^ = cxd (chain-shift model) and the model 
of three-dimensional finite-size defects (^ = 30 A). 
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tion of defects with a given concentration p. Assuming 
different dimensionalities and correlation lengths of the 
triclinic islands, the field distribution D(B^) and the cor- 
responding muon relaxation functions were calculated, 
Pit) = 3 + |E'=i/j/^(^Mj)cos(7^B^,,t)dS^,,. In 
Fig- ESI field distributions at the two muon stopping sites 
are compared for the model of one-dimensional defects 
with infinite ^ and the model of three-dimensional defects 
with the size ^ = 30 A, both calculated for p = 0.35. The 
former yields a few sharp peaks and the latter a broader 
field distribution. The two-peak field-distribution at the 
^^Na site obtained for the second model is also shown. 
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